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SOME EFFECTS OF TIP FINS ON WING FLUTTER CHARACTERISTICS 


By Robi>rt C. Goetz and Robert V. Doggett, Jr. 

Langley Research Center 

SUMMARY 

A wind-tunnel investigation has been conducted over the Mach number range from 
about 0.6 to 1.2 to determine the effects of large tip fins on the flutter characteristics of 
a swept wing. The basic wing configuration had an aspect ratio of 0.95, leading -edge 
sweep of 40*^, and trailing -edge sweep of 21*^. Two of these configurations were modified 
with tip fins of 60*^ dihedral and had effective aspect ratios of 1.5 and 2.2. In general, 
the results indicate that the addition of tip fins reduces the flutter speed, with the larger 
fin having the greater effect. 

Comparison of the experimental flutter speeds at Mach numbers between 0.60 and 
0.90 with calculated values obtained by using doublet -lattice unsteady aerodynamic theory 
was good. Analytical results where structural and aerodynamic effects of the tip fins 
were isolated indicated that the reduction in flutter speed produced by the addition of the 
fins was caused by both effects, with the structural effect being the more pronounced. 

INTRODUCTION 

For the vehicle designs emerging for space transportation systems, it is evident 
that interactive effects from both a steady and unsteady aerod}mamic viewpoint are 
becoming more significant. These vehicle designs are combinations of multiple bodies 
and lifting surfaces; the latter are dictated by requireme»'«’c for the maneuverability and 
the stability and control during the ascent, reentry, glidu.,', and landing phases of flight. 
Lifting-surface concepts that appear attractive for satisfying these requirements may 
have the feature of noi^lanar intersecting surface segments. One such concept is a 
swept wing with large tip fins. 

Wings with large tip fins are a concern to the aeroelastician since their structural 
design might be significantly influenced by flutter clearance requirements. For example, 
if the mass of the tip fin is not negligible in relation to that of the wing alone, the flutter 
speed is modified by its presence. Whereas a change in mass distribution is known to 
have a large effect on the flutter speed, the associated change in aerodynamic forces may 
also be a significant factor. 
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Concern over designs where aerodynamic interference may be important to flutter 
is not new. When T-tail designs began to emerge in the early 1950' s, it was quickly rec- 
ognized that a better understanding of aerodynamic interference was needed to determine 
the flutter characteristics of these aerodynamically complex designs. This need resulted 
in the development of unsteady aerodynamic theories ^plicable to T -tails such as those 
presented in references 1 and 2. In addition to analytical studies some experimental 
T-tail flutter results are available. (For example, see refs. 3 and 4.) Another area that 
has received attention is tandem -mounted wing — horizontal-tail configurations. Several 
wind-tunnel -model flutter studies have been conducted on wing -tail configurations, and 
some results are presented in references 5 and 6. Unsteady aerodynamic theories have 
been developed for application to wing -tail configurations. A kernel -function procedure 
is presented in reference 7, and the doublet -lattice method is described in reference 8. 
Some experimental flutter results are compared in reference 6 with anal 3 dical results 
obtained by using both kernel -function and doublet -lattice unsteady aerodynamic methods. 
Additional comparisons between experimental and analytical results are presented in 
reference 9. However, a search of the available literature indicated that there are no 
published flutter results in the range from high subsonic to low supersonic speeds for 
wings with large tip fins. 

Accordingly, the present investigation was undertaken to determine the importance 
of the aerodynamic interference effects on the flutter characteristics of a swept wing hav- 
ing a tip fin and to compare experimental and analytical results for subsonic flow obtained 
from application of the doublet -lattice method. 

The experimental program was conducted by using simple semispan models in the 
Langley 26-inch transonic tunnel over the Mach number range from about 0.6 to 1.2. 

Three model configurations were tested. The basic wing -alone configuration had an 
aspect ratio of 0.95, leading-edge sweep of 40°, and trailing -edge sweep of 21°. The 
other two model configurations had the same leading- and trailing -edge sweep angles, but 
tip fins with 60° of dihedral (with no toe-in) were added so that the resulting models had 
effective aspect ratios of 1.5 and 2.2. All model configurations had 10 -percent -thick 
NACA 64 -series airfoil sections and were tested at zero angle of attack. A coupled -mode 
flutter analysis based on doublet -lattice unsteady aerodynamics was performed for Mach 
numbers of 0.60, 0.80, and 0.90 by using calculated mode shapes and measured frequen- 
cies; the results of this analysis are compared with measured results. 

SYMBOLS 

br reference semichord, m 

f frequency, Hz 
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flutter frequency, Hz 

frequency of second natural mode, Hz 

structural damping coefficient 

Mach number 

mass, kg 

dynamic pressure, -^pV^, kN/m2 = kPa 
velocity, meters/second 
flutter-speed index parameter, - — — 

reference volume, m^ 
mass -ratio parameter, ^ 
density, kg/m3 

reference circular frequency, 2 ?rf 2 , rad/sec 


Subscripts: 
c calculated 

e experimental 


MODELS 

Description 

The basic model configuration (wing alone) used in this investigation was a semi- 
span, aspect-ratio-0.95 wing with no dihedral and having leading -edge sweep of 40° and 
trailing -edge sweep of 21°. Two other models, consisting of tho basic configuration with 
two different lengths of wing -tip extensions added at a dihedral angle of 60°, were tested 
in an effort to determine the effects of a tip fin on the flutter characteristics of the basic 
wing. The fins had the same leading- and trailing -edge sweep as the basic wing and were 



oriented at zero angle of attack with respect to the free -stream flow. The two fin mod- 
els, referred to hereafter as the full -fin and half -fin models, had effective aspect ratios 
of 2.2 and 1.5, respectively. The effective aspect ratio is defined as the square of the 
semispan normalized by the total wing area after the fin has been rotated into the plane 
of the wing. Both the wing and the fins had 10 -percent -thick NACA 64 -series airfoil 
sections parallel to the flow. The fins were relatively large compared with the wing; 
for the full -fin model the fin planform area was about half that of the basic wing, and for 
the half -fin model the fin planform area was slightly less than one -third that of the basic 
wing. 

Figure 1 is a photograph illustrating the three model configurations. A sketch 
giving the geometric properties of the models is presented in figure 2. 

The models were of simple sandwich construction, consisting of a 0.2235-cm-thick 
aluminum -alloy core to which balsa wood was bonded with the grain oriented perpendic- 
ular to the core. The balsa wood was machined to give the desired NACA 64 AO 10 airfoil 
shape. The aluminum -alloy core of the model was also used to support the model in a 
cantilever fashion along the forward two-thirds of the model root. Resistance -type 
strain gages were mounted on the model core to measure model dynamic response. 

Physical Properties and Vibration Characteristics 

The total measured mass, natural frequencies f, and structural damping g of 
the models are given in table I. The natural frequencies, node lines, and still -air damp- 
ing coefficients of the first four natural structural modes were determined for all three 
models. The natural frequency for the full -fin -model fifth mode was also determined. 
The calculated and measured frequencies and nodal patterns are shown in figure 3. 

These experimental data were obtained by exciting the models at resonance with an inter- 
rupter air -jet shaker using the Ig sand technique. For these tests the models were 
attached to a massive backstop. The mounting arrangement was essentially the same as 
that used for mounting the models in the wind tunnel. The damping data were obtained 
from the decaying oscillation that resulted when the pulsating air jet from the shaker was 
abruptly cut off. The calculated modal data were obtained 1 y using the NASA Structural 
Analysis (NASTRAN) Computer Program. NASTRAN is described in detail in refer- 
ences 10 and 11. Quadrilateral structural finite elements (NASTRAN QUAD2) were used 
to model the structure for the calculations. Thirty-six elements were used for the wing 
portion of all three models; 24 elements were used for the half -tip-fin portion; and 
60 elements were used for the full -tip -fin portion. The arrangement of the quadrilateral 
elements is shown in figure 4. A comparison of the calculated and measured node lines 
(see fig. 3) indicates fairly good agreement. Although all the frequencies agree reason- 
ably well, ordinarily the best agreement would have been expected for the basic wing 
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model since this is the simplest structural configuration. However, better agreement 
was obtained for the two models with tip fins which are structurally more complex than 
the wing model. The best frequency agreement was obtained for the full -fin model. 
Although no specific explanation is available, it should be pointed out that some uncer- 
tainty was introduced into the analytical model since it was necessary to estimate the 
stiffness properties of the balsa wood and of the adhesive used to attach the balsa -wood 
covering to the aluminum -alloy core. The stiffness properties of the aluminum were 
known, but it was necessary to estimate the properties of the balsa wood and the adhe- 
sive by use of typical values of similar material. The estimated values of Young's mod- 
ulus used for the balsa wood and for the adhesive were 2.83 GPa and 3.45 GPa, respec- 
tively. The adhesive was assumed to be 0.056 cm thick. 

FLUTTER EXPERIMENTS 
Wind Tunnel 

The investigation was conducted in the Langley 26 -inch transonic tunnel over the 
Mach number range from about 0.6 to 1.2. This transonic blowdown tunnel has an octag- 
onal test section which measures 66 cm across the flats and has a slot in each corner. 
The tunnel is capable of qieration at stagnation pressures up to about 517 kPa. During 
operation of the tunnel, the area of the orifice (second minimum) is preset at a given 
value; as the stagnation pressure, and thus the density, is increased, the test -section 
Mach number increases until the orifice becomes choked. Thereafter, as the stagnation 
pressure is further increased, the Mach number remains approximately constant. How- 
ever, the area of the orifice may also be varied during a test run as the stagnation pres- 
sure is increased, or held constant, so that various operating paths of Mach number and 
density may be followed. Both methods of operation were used in the present 
investigation. 


Test Procedure 

The models tested were mounted at midheight on a 7.6-cm-diameter fuselage 
sting (see fig. 5) which extended forward into the low -speed region of the tunnel. This 
arrangement prevented the formation of shock waves from the sting nose which might 
reflect from the tunnel walls onto the models. The sting provided a rigid mount for the 
models since the mass of the support system was very large compared with the mass of 
a model. The fundamental frequency of the support system was about 6 Hz, well below 
the first natural frequency of any of the models. 

An optical system displayed an image of the model on a ground-glass screen dur- 
ing the wind-tunnel tests. When flutter was observed visually, the airflow was quickly 
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stopped in an effort to save the model from being damaged so that it could be utilized in 
subsequent tests. Before each test run, with the model mounted in the tunnel, frequen- 
cies and damping coefficients for the first two natural modes were measured to be cer- 
tain that the model had experienced no structural damage. The tunnel stagnation pres- 
sure, stagnation temperature, test -section static pressure, and model strain -gage 
signals were continuously recorded on a direct readout recorder. Visual records of the 
model motion were obtained through the use of two high-speed motion-picture cameras. 

FLUTTER ANALYSIS 

Subsonic flutter calculations were made for all three models over the Mach number 
range from 0.60 to 0.90. The flutter equations in matrix notation were expressed in 
terms of generalized modal coordinates, and the traditional V -g method of solution was 
used. The calculated natural mode shapes were used in conjunction with the correspond- 
ing measured natural frequencies. 

Flutter characteristics for all three models were determined by using doublet - 
lattice unsteady aerodynamic theory. The doublet -lattice unsteady aerodynamic forces 
were determined by using the method described in references 12 and 13. The doublet - 
lattice method requires that the lifting surface be subdivided into trapezoidal boxes 
arranged in streamwise columns. The arrangement of boxes used is shown in figure 6. 
The same arrangement was used for the wing portion of all three models. There were 
50 boxes on the wing, 35 boxes on the half fin, and 50 boxes on the full fin. The two fins 
shown in the figure have been rotated down into the plane of the paper so that the view 
is normal to these surfaces. For the doublet -lattice method, a line of acceleration 
potential doublets is placed at the one -quarter -chord station of each box. An aerody- 
namic influence coefficient matrix is generated which relates the force on the boxes to 
the downwash on the boxes. The force acts at the one -quarter -chord point, and the 
downwash point where the geometrical boundary condition of tangential flow is satisfied 
is at the three -quarter -chord station. Both the force and downwash points are located 
at the box midspan station. The aerodynamic influence coefficient is used in conjunction 
with the mode shapes to determine a generalized aerodynamic force matrix. 

Subsonic lifting-surface theory (kernel function) was also used to calculate the 
flutter characteristics of the basic wing model at a Mach number of 0.80. The technique 
used to generate the kernel-function unsteady aerodynamic forces was based on that 
described in reference 14. For the kernel -function method, a linear integral equation 
which relates the pressure distribution to the downwash at selected control points on the 
lifting surface is solved numerically for the unsteady pressure distribution. The down- 
wash at the control points is a function of the structural mode shapes used in the analy- 
sis. The kernel -function method generates directly a generalized unsteady aerodynamic 
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force matrix. Thirty-six control points were used in the kernel -function analysis. The 
locations of these points are indicated by the circular symbols in figure 6. 

Since the aerodynamic models for the doublet -lattice and kernel -function methods 
were different from the structural model (compare figs. 4 and 6), it was necessary to 
Interpolate the modal deformations determined at the structural grid points to the modal 
displacements and streamwise slopes required at the aerodynamic points. This trans- 
formation was accomplished by using the surface spline function described in refer- 
ence IS. Since the surface spline function is based on the small deflection equation of 
an infinite plate, it is very good for interpolating structural mode shapes. In interpo- 
lating the modes for the two models with tip fins, two spline functions were used. One 
spline function was generated by using the modal displacements on the wing portion, and 
the second was generated by using the modal displacements on the fin portion. The 
modal deformations along the juncture between the wing and fin portions were included 
in the generation of both splines to assure continuity in the interpolation. 

The V-g method of flutter solution requires many solutions of a complex eigenvalue 
problem with the reduced frequency as a parameter and is relatively expensive in terms 
of computer costs. Since one of the most costly processes is the determination of the 
generalized unsteady aerodynamic forces, it has become more or less standard practice 
to calculate the generalized aerodynamic forces for a relatively small number of values 
of the reduced frequency parameter and to interpolate these forces for their values at a 
larger number of reduced frequencies. Interpolated aerodynamic forces were used to 
obtain the calculated flutter results presented in this paper. A natural cubic spline was 
used to perform the interpolation. A discussion of the use of cubic spline functions for 
aerodynamic force interpolation is presented in reference 16. For the doublet -lattice 
calculations, the aerodynamic forces were calculated at 6 values of reduced frequency 
and interpolated to 50 values of reduced frequency. For the kernel function, there were 
20 calculated values and 400 interpolated values. 

RESULTS AND DISCUSSION 

The basic experimental flutter results are presented in table n and figure 7. The 
data presented in the figure include the variations with Mach number of the mass -ratio 
parameter 4, of the flutter -frequency ratio ff/f2> ^nd of the flutter-speed index param- 
eter Vj. The flutter-speed-index-parameter curves represent stability boundaries, 
with the stable region (no flutter) below the curve. This parameter depends on the phys- 
ical properties of the model, in particular the stiffness, and the atmosphere in which it 
operates. When plotted as the ordinate against Mach number, curves of constant dynamic 
pressure are lines parallel to the Mach number abscissa. The mass-ratio parameter 
4 is defined as the ratio of the mass of the model to the mass of a representative 
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surrounding volume of test medium. The volume used in this study is that contained in 
the conical frustum having base diameters equal to the respective model root and tip 
chords and a height equal to the sum of the individual spans of the wing and tip fin por- 
tions where the span of the fin is measured in the plane of the fin. These volumes for 
the wing, half -fin, and full -fin models were 2061.9 cm3, 2504.5 cm3, and 2740.6 cm3, 
respectively. The second natural frequency was used as the reference frequency. 
The same reference length br was used for all three models. This length 6.26 cm 
was the semichord at the three-quarter-span station of the wing portion of the models. 

No unusual trends are shown by the data presented in figure 7 for the half -fin and 
full-fin models. The flutter boundary is similar to that usually observed, namely, a 
more or less constant value at subsonic speeds followed by a decrease in flutter -speed 
index at transonic Mach numbers with a sharp increase in flutter -speed index as the 
flow becomes supersonic. No experimental flutter results were obtained for the basic 
wing model. This n?.odel did not flutte v within the available (derating range of the wind 
tunnel. At transonic speeds this model was tested at dynamic pressure levels of about 
140 kPa which is considerably above the flutter boundary determined for the other two 
models. 

In terms of the flutter -speed index parameter, a comparison of the flutter bound- 
ary for the half -fin model (fig. 7(a)) with that of the full-fin model (fig. 7(b)) would lead 
the reader to believe that the full-fin configuration is less susceptible to flutter. This 
is not the true situation, however, since the full-fin model fluttered at lower dynamic 
pressures than did the half -fin model over the Mach number range of the study. Nor- 
mally the flutter-speed index parameter Vj is successful in correlating flutter data 
between similar configurations; however, this was not true for the wing -tip configura- 
tions for this study. The variation of flutter dynamic pressure with Mach number for 
both of these models is presented in figure 8. Both sets of data show the same trend, 
namely, a more or less constant value at subsonic speeds, a dramatic decrease in the 
transonic range, and an increase in value at supersonic speeds. For both models, the 
transonic minimum dynamic pressure was about 45 percent of the value at Mach 0.6. 
Although no flutter data were obtained for the basic wing model, it is known that its flut- 
ter boundary is considerably above those of the models with fins. As mentioned previ- 
ously, the wing model was tested at transonic Mach numbers to a dynamic pressure of 
140 kPa. Therefore, it can be concluded that the addition of the fins to the basic wing 
has a significant adverse effect on the flutter characteristics. Further, the larger the 
fin, the greater is the decrease in flutter speed. 

Experimental and anal 3 rtical flutter boundaries for the two models with tip fins are 
compared in figure 9. The data are presented as the variation with Mad. number of the 
flutter-frequency ratio and the flutter-speed index parameter. The experim<;ntal curves 
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shown are the faired curves from figure 7. Additional analytical flutter results, includ- 
ing basic wing model data, are presented in table ni. The calculated results in the fig- 
ure were obtained by using a four -mode analysis with doublet -’ettice unsteady aerody- 
namic forces and the densities corresponding to the expevi:. .«ittm' alues. One calculation 
was made at M » 0.80 (see table ni(c)) for the full -fin ..lOdel by l.- ig five natural modes 
in the analysis. Since this result was almost the same as the four-mode result, it was 
assumed that the use of four modes was sufficient to insure convergence of the flutter 
solution. Also, some calculated flutter results vere obtained by using densities differ- 
ent from the e]Q)erimental values to assess the sensitivity of the calculated flutter speeds 
to mass -ratio effects. Although the calculated flutter speeds did indicate a sensitivity of 
flutter speed to changes in density, no large effect was observed over the range of the 
study. In general, the calculated flutter results predict higher flutter frequencies and 
speeds than were found experimentally. However, the agreement between the two sets of 
data is considered to be good and shows that the doublet -lattice unsteady aerodynamic 
theory can be used with some confidence in predicting the flutter speeds of wings with 
large tip fins. 

As has been observed, there was a marked reduction in the flutter speed when a 
tip fin was added to the basic wing model. This effect could be caused by both structural 
and aerodynamic effects. An analytical study was made to separate the aerodynamic and 
structural effects, and the results of this study are shown in figure 10. The data are 
presented in the form of the variation of flutter dynamic pressure with Mach number for 
the wing model, the full -fin model with no aerodynamic forces acting on the fin portion, 
and the full -fin model with aerodynamic forces acting on the fin. Additional calculated 
flutter results for all three cases are presented in table HI. All results presented in the 
figure were obtained by using the doublet -lattice unsteady aerodynamic theory with the 
exception of a M » 0.80 result obtained by using kernel -function aerodynamics for the 
wing model. It is interesting to note that the kernel -function and doublet -lattice results 
are almost the same. The calculated full-fin model results are the same as those pre- 
sented in figure 9(b). Flutter characteristics for the full-fin model without aerodynamic 
forces on the fin (hereafter referred to as full fin without aerodynamics) were calculated 
by using a five -mode analysis and the densities corresponding to the full-fin model exper- 
iment. A check calculation at M = 0.80 using only four modes for this configuration 
gave the same result as was obtained by using five modes. (See table 111(d).) Conse- 
quently, it can be assumed that the flutter solutions using five mode.-i for the full fin with- 
out aerodynamics case have converged. The basic-wing -model results presented in the 
figure were obtained by using four modes and an air density of 2.701 kg/m^. For both 
the basic wing model and the full fin without aerod)mamics case, additional calculations 
were made with density as a variable to determine whether any significant mass -ratio 
effects were present. (See tables Ill(a) and m(d).) No significant mass-ratio effects 
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were noted over the density range covered. An examination of the data in figure 10 
(compare upper and lower curves) indicates that there is a very large decrease in the 
flutter dynamic pressure with the addition of the full fin to the basic wing model. For 
example, at M = 0.60 this reduction is 77 percent and includes both aerodynamic and 
structural effects. By removing the aerodynamic forces from the full fin, it is possible 
to determine what proportion of the reduction in flutter dynamic pressure is due to the 
structural changes which occur when the fin is added to the wing. The flutter boundary 
without fin aerodynamic forces is shown as the middle curve in the figure. This flutter 
boundary is substantially lower than the basic wing boundary. At M = 0.60 the flutter 
dynamic pressure is reduced by 63 percent. The results in figure 10 incJ ;ate that for 
the configuration studied in this investigation, the greatest proportion of the reduction in 
flutter dynamic pressure was due to the structural and inertial changes p^'oduced by the 
addition of the tip fin. However, although the aerodynamic effects were t-maller than the 
structural effects, they were not insignificant. The aerodynamic effect was responsible 
for about 20 percent of the total decrease in flutter dynamic pressure. It should be 
noted that the relative proportion of structural and aerodynamic effects of tip fins on 
wing flutier characteristics would differ for other configurations. 

CONCLUDING REMARKS 

The effects of the addition of two different tip fins >asic wing have been deter- 
mined experimentally over the Mach munber range frot jt 0.6 to 1.2. The basic 
wing had an aspect ratio of 0.9S, a leading -edge sweep oi tO^ and a trailing -edge sweep 
of 21^. Both the small and large tip fins had a dihedral of 60^ and, in terms of surface 
area, were one -third and one -half as large as the basic wing, respectively. The re 
indicate that the addition of the tip fins reduces the flutter speed, with the larger fin 
having the greater effect. No unusual behavior was noted in the variation of the flutter 
speed with Mach number for the models with tip fins. 

Comparison of the experimental flutter boundaries at subsonic speeds with calcu- 
lated results obtained by using doublet -lattice unsteady aerodynamic theory was good. 
The results of an analytical study where structural and aerodynamic effects of the tip 
fins were isolated indicated that the reduction in flutter speed produced by the addition of 
the fins was caused by a combination of structural and aerodynamic effects, with the 
structural effect being the more pronounced. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 12, *974. 
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Figure 10.- Calculated effect of tip fin on wing flutter d 3 rnamic pressure. 


